Assessing transmission risks and control strategy for monkeypox as an emerging zoonosis in a metropolitan area.

To model the spread of monkeypox (MPX) in a metropolitan area for assessing the risk of possible outbreaks, and identifying essential public health measures to contain the virus spread. The animal reservoir is the key element in the modeling of zoonotic disease. Using a One Health approach, we model the spread of the MPX virus in humans considering potential animal hosts such as rodents (e.g., rats, mice, squirrels, chipmunks, etc.) and emphasize their role and transmission of the virus in a high-risk group, including gay and bisexual men-who-have-sex-with-men (gbMSM). From model and sensitivity analysis, we identify key public health factors and present scenarios under different transmission assumptions. We find that the MPX virus may spill over from gbMSM high-risk groups to broader populations if the efficiency of transmission increases in the higher-risk group. However, the risk of outbreak can be greatly reduced if at least 65% of symptomatic cases can be isolated and their contacts traced and quarantined. In addition, infections in an animal reservoir will exacerbate MPX transmission risk in the human population. Regions or communities with a higher proportion of gbMSM individuals need greater public health attention. Tracing and quarantine (or "effective quarantine" by postexposure vaccination) of contacts with MPX cases in high-risk groups would have a significant effect on controlling the spreading. Also, monitoring for animal infections would be prudent.

Modeling vaccination and control strategies for outbreaks of monkeypox at gatherings.

Background: The monkeypox outbreak in non-endemic countries in recent months has led the World Health Organization (WHO) to declare a public health emergency of international concern (PHEIC). It is thought that festivals, parties, and other gatherings may have contributed to the outbreak. Methods: We considered a hypothetical metropolitan city and modeled the transmission of the monkeypox virus in humans in a high-risk group (HRG) and a low-risk group (LRG) using a Susceptible-Exposed-Infectious-Recovered (SEIR) model and incorporated gathering events. Model simulations assessed how the vaccination strategies combined with other public health measures can contribute to mitigating or halting outbreaks from mass gathering events. Results: The risk of a monkeypox outbreak was high when mass gathering events occurred in the absence of public health control measures. However, the outbreaks were controlled by isolating cases and vaccinating their close contacts. Furthermore, contact tracing, vaccinating, and isolating close contacts, if they can be implemented, were more effective for the containment of monkeypox transmission during summer gatherings than a broad vaccination campaign among HRG, when accounting for the low vaccination coverage in the overall population, and the time needed for the development of the immune responses. Reducing the number of attendees and effective contacts during the gathering could also prevent a burgeoning outbreak, as could restricting attendance through vaccination requirements. Conclusion: Monkeypox outbreaks following mass gatherings can be made less likely with some restrictions on either the number and density of attendees in the gathering or vaccination requirements. The ring vaccination strategy inoculating close contacts of confirmed cases may not be enough to prevent potential outbreaks; however, mass gatherings can be rendered less risky if that strategy is combined with public health measures, including identifying and isolating cases and contact tracing. Compliance with the community and promotion of awareness are also indispensable to containing the outbreak.

Community structured model for vaccine strategies to control COVID19 spread: A mathematical study.

Initial efforts to mitigate the COVID-19 pandemic have relied heavily on non-pharmaceutical interventions (NPIs), including physical distancing, hand hygiene, and mask-wearing. However, an effective vaccine is essential to containing the spread of the virus. We developed a compartmental model to examine different vaccine strategies for controlling the spread of COVID-19. Our framework accounts for testing rates, test-turnaround times, and vaccination waning immunity. Using reported case data from the city of Toronto, Canada between Mar-Dec, 2020 we defined epidemic phases of infection using contact rates as well as the probability of transmission upon contact. We investigated the impact of vaccine distribution by comparing different permutations of waning immunity, vaccine coverage and efficacy throughout various stages of NPI's relaxation in terms of cases and deaths. The basic reproduction number is also studied. We observed that widespread vaccine coverage substantially reduced the number of cases and deaths. Under phases with high transmission, an early or late reopening will result in new resurgence of the infection, even with the highest coverage. On the other hand, under phases with lower transmission, 60% of coverage is enough to prevent new infections. Our analysis of R0 showed that the basic reproduction number is reduced by decreasing the tests turnaround time and transmission in the household. While we found that household transmission can decrease following the introduction of a vaccine, public health efforts to reduce test turnaround times remain important for virus containment.

Mathematical modelling of vaccination rollout and NPIs lifting on COVID-19 transmission with VOC: a case study in Toronto, Canada.

BACKGROUND: Since December 2020, public health agencies have implemented a variety of vaccination strategies to curb the spread of SARS-CoV-2, along with pre-existing Nonpharmaceutical Interventions (NPIs). Initial strategies focused on vaccinating the elderly to prevent hospitalizations and deaths, but with vaccines becoming available to the broader population, it became important to determine the optimal strategy to enable the safe lifting of NPIs while avoiding virus resurgence. METHODS: We extended the classic deterministic SIR compartmental disease-transmission model to simulate the lifting of NPIs under different vaccine rollout scenarios. Using case and vaccination data from Toronto, Canada between December 28, 2020, and May 19, 2021, we estimated transmission throughout past stages of NPI escalation/relaxation to compare the impact of lifting NPIs on different dates on cases, hospitalizations, and deaths, given varying degrees of vaccine coverages by 20-year age groups, accounting for waning immunity. RESULTS: We found that, once coverage among the elderly is high enough (80% with at least one dose), the main age groups to target are 20-39 and 40-59 years, wherein first-dose coverage of at least 70% by mid-June 2021 is needed to minimize the possibility of resurgence if NPIs are to be lifted in the summer. While a resurgence was observed for every scenario of NPI lifting, we also found that under an optimistic vaccination coverage (70% coverage by mid-June, along with postponing reopening from August 2021 to September 2021) can reduce case counts and severe outcomes by roughly 57% by December 31, 2021. CONCLUSIONS: Our results suggest that focusing the vaccination strategy on the working-age population can curb the spread of SARS-CoV-2. However, even with high vaccination coverage in adults, increasing contacts and easing protective personal behaviours is not advisable since a resurgence is expected to occur, especially with an earlier reopening.

Efficacy of a "stay-at-home" policy on SARS-CoV-2 transmission in Toronto, Canada: a mathematical modelling study.

BACKGROUND: Globally, nonpharmaceutical interventions for COVID-19, including stay-at-home policies, limitations on gatherings and closure of public spaces, are being lifted. We explored the effect of lifting a stay-at-home policy on virus resurgence under different conditions. METHODS: Using confirmed case data from Toronto, Canada, between Feb. 24 and June 24, 2020, we ran a compartmental model with household structure to simulate the impact of the stay-at-home policy considering different levels of compliance. We estimated threshold values for the maximum number of contacts, probability of transmission and testing rates required for the safe reopening of the community. RESULTS: After the implementation of the stay-at-home policy, the contact rate outside the household fell by 39% (from 11.58 daily contacts to 7.11). The effective reproductive number decreased from 3.56 (95% confidence interval [CI] 3.02-4.14) on Mar. 12 to 0.84 (95% CI 0.79-0.89) on May 6. Strong adherence to stay-at-home policies appeared to prevent SARS-CoV-2 resurgence, but extending the duration of stay-at-home policies beyond 2 months had little added effect on cumulative cases (25 958 for 65 days of a stay-at-home policy and 23 461 for 95 days, by July 2, 2020) and deaths (1404 for 65 days and 1353 for 95 days). To avoid a resurgence, the average number of contacts per person per day should be kept below 9, with strict nonpharmaceutical interventions in place. INTERPRETATION: Our study demonstrates that the stay-at-home policy implemented in Toronto in March 2020 had a substantial impact on mitigating the spread of SARS-CoV-2. In the context of the early pandemic, before the emergence of variants of concern, reopening schools and workplaces was possible only with other nonpharmaceutical interventions in place.

Introduction to Focus Issue: Dynamical disease: A translational approach.

The concept of Dynamical Diseases provides a framework to understand physiological control systems in pathological states due to their operating in an abnormal range of control parameters: this allows for the possibility of a return to normal condition by a redress of the values of the governing parameters. The analogy with bifurcations in dynamical systems opens the possibility of mathematically modeling clinical conditions and investigating possible parameter changes that lead to avoidance of their pathological states. Since its introduction, this concept has been applied to a number of physiological systems, most notably cardiac, hematological, and neurological. A quarter century after the inaugural meeting on dynamical diseases held in Mont Tremblant, Québec [Bélair et al., Dynamical Diseases: Mathematical Analysis of Human Illness (American Institute of Physics, Woodbury, NY, 1995)], this Focus Issue offers an opportunity to reflect on the evolution of the field in traditional areas as well as contemporary data-based methods.

Fangcang shelter hospitals during the COVID-19 epidemic, Wuhan, China.

OBJECTIVE: To design models of the spread of coronavirus disease-2019 (COVID-19) in Wuhan and the effect of Fangcang shelter hospitals (rapidly-built temporary hospitals) on the control of the epidemic. METHODS: We used data on daily reported confirmed cases of COVID-19, recovered cases and deaths from the official website of the Wuhan Municipal Health Commission to build compartmental models for three phases of the COVID-19 epidemic. We incorporated the hospital-bed capacity of both designated and Fangcang shelter hospitals. We used the models to assess the success of the strategy adopted in Wuhan to control the COVID-19 epidemic. FINDINGS: Based on the 13 348 Fangcang shelter hospitals beds used in practice, our models show that if the Fangcang shelter hospitals had been opened on 6 February (a day after their actual opening), the total number of COVID-19 cases would have reached 7 413 798 (instead of 50 844) with 1 396 017 deaths (instead of 5003), and the epidemic would have lasted for 179 days (instead of 71). CONCLUSION: While the designated hospitals saved lives of patients with severe COVID-19, it was the increased hospital-bed capacity of the large number of Fangcang shelter hospitals that helped slow and eventually stop the COVID-19 epidemic in Wuhan. Given the current global pandemic of COVID-19, our study suggests that increasing hospital-bed capacity, especially through temporary hospitals such as Fangcang shelter hospitals, to isolate groups of people with mild symptoms within an affected region could help curb and eventually stop COVID-19 outbreaks in communities where effective household isolation is not possible.

On the Effect of Age-Dependent Mortality on the Stability of a System of Delay-Differential Equations Modeling Erythropoiesis.

We present an age-structured model for erythropoiesis in which the mortality of mature cells is described empirically by a physiologically realistic probability distribution of survival times. Under some assumptions, the model can be transformed into a system of delay differential equations with both constant and distributed delays. The stability of the equilibrium of this system and possible Hopf bifurcations are described for a number of probability distributions. Physiological motivation and interpretation of our results are provided.

Transit and lifespan in neutrophil production: implications for drug intervention.

A comparison of the transit compartment ordinary differential equation modelling approach to distributed and discrete delay differential equation models is studied by focusing on Quartino's extension to the Friberg transit compartment model of myelosuppression, widely relied upon in the pharmaceutical sciences to predict the neutrophil response after chemotherapy, and on a QSP delay differential equation model of granulopoiesis. An extension to the Quartino model is provided by considering a general number of transit compartments and introducing an extra parameter that allows for the decoupling of the maturation time from the production rate of cells. An overview of the well established linear chain technique, used to reformulate transit compartment models with constant transit rates as distributed delay differential equations (DDEs), is then given. A state-dependent time rescaling of the Quartino model is performed to apply the linear chain technique and rewrite the Quartino model as a distributed DDE, yielding a discrete DDE model in a certain parameter limit. Next, stability and bifurcation analyses are undertaken in an effort to situate such studies in a mathematical pharmacology context. We show that both the original Friberg and the Quartino extension models incorrectly define the mean maturation time, essentially treating the proliferative pool as an additional maturation compartment. This misspecification can have far reaching consequences on the development of future models of myelosuppression in PK/PD.

Normal and pathological dynamics of platelets in humans.

We develop a mathematical model of platelet, megakaryocyte, and thrombopoietin dynamics in humans. We show that there is a single stationary solution that can undergo a Hopf bifurcation, and use this information to investigate both normal and pathological platelet production, specifically cyclic thrombocytopenia. Carefully estimating model parameters from laboratory and clinical data, we then argue that a subset of parameters are involved in the genesis of cyclic thrombocytopenia based on clinical information. We provide model fits to the existing data for both platelet counts and thrombopoietin levels by changing four parameters that have physiological correlates. Our results indicate that the primary change in cyclic thrombocytopenia is an interference with, or destruction of, the thrombopoietin receptor with secondary changes in other processes, including immune-mediated destruction of platelets and megakaryocyte deficiency and failure in platelet production. This study contributes to the understanding of the origin of cyclic thrombocytopenia as well as extending the modeling of thrombopoiesis.

Walking on a Vertically Oscillating Treadmill: Phase Synchronization and Gait Kinematics.

Sensory motor synchronization can be used to alter gait behavior. This type of therapy may be useful in a rehabilitative setting, though several questions remain regarding the most effective way to promote and sustain synchronization. The purpose of this study was to describe a new technique for using synchronization to influence a person's gait and to compare walking behavior under this paradigm with that of side by side walking. Thirty one subjects walked on a motorized treadmill that was placed on a platform that oscillated vertically at various frequencies and amplitudes. Synchronization with the platform and stride kinematics were recorded during these walking trials and compared with previously reported data from side by side walking. The results indicated that vertical oscillation of the treadmill surface at frequencies that matched subjects preferred stride or step frequency resulted in greater unintentional synchronization when compared with side by side walking data (up to 78.6±8.3% of the trial vs 59.2±17.4%). While intermittent phase locking was observed in all cases, periods of synchronization occurred more frequently and lasted longer while walking on the oscillating treadmill (mean length of periods of phase locking 11.85 steps vs 5.18 steps). Further, stride length, height and duration were altered by changing the frequency of treadmill oscillation. These results suggest that synchronization to a haptic signal may hold implications for use in a clinical setting.

Modeling the spread and control of dengue with limited public health resources.

A deterministic model for the transmission dynamics of dengue fever is formulated to study, with a nonlinear recovery rate, the impact of available resources of the health system on the spread and control of the disease. Model results indicate the existence of multiple endemic equilibria, as well as coexistence of an endemic equilibrium with a periodic solution. Additionally, our model exhibits the phenomenon of backward bifurcation. The results of this study could be helpful for public health authorities in their planning of a proper resource allocation for the control of dengue transmission.

Neutrophil dynamics during concurrent chemotherapy and G-CSF administration: Mathematical modelling guides dose optimisation to minimise neutropenia.

The choice of chemotherapy regimens is often constrained by the patient's tolerance to the side effects of chemotherapeutic agents. This dose-limiting issue is a major concern in dose regimen design, which is typically focused on maximising drug benefits. Chemotherapy-induced neutropenia is one of the most prevalent toxic effects patients experience and frequently threatens the efficient use of chemotherapy. In response, granulocyte colony-stimulating factor (G-CSF) is co-administered during chemotherapy to stimulate neutrophil production, increase neutrophil counts, and hopefully avoid neutropenia. Its clinical use is, however, largely dictated by trial and error processes. Based on up-to-date knowledge and rational considerations, we develop a physiologically realistic model to mathematically characterise the neutrophil production in the bone marrow which we then integrate with pharmacokinetic and pharmacodynamic (PKPD) models of a chemotherapeutic agent and an exogenous form of G-CSF (recombinant human G-CSF, or rhG-CSF). In this work, model parameters represent the average values for a general patient and are extracted from the literature or estimated from available data. The dose effect predicted by the model is confirmed through previously published data. Using our model, we were able to determine clinically relevant dosing regimens that advantageously reduce the number of rhG-CSF administrations compared to original studies while significantly improving the neutropenia status. More particularly, we determine that it could be beneficial to delay the first administration of rhG-CSF to day seven post-chemotherapy and reduce the number of administrations from ten to three or four for a patient undergoing 14-day periodic chemotherapy.

Threshold dynamics in an SEIRS model with latency and temporary immunity.

A disease transmission model of SEIRS type with distributed delays in latent and temporary immune periods is discussed. With general/particular probability distributions in both of these periods, we address the threshold property of the basic reproduction number R0 and the dynamical properties of the disease-free/endemic equilibrium points present in the model. More specifically, we 1. show the dependence of R0 on the probability distribution in the latent period and the independence of R0 from the distribution of the temporary immunity, 2. prove that the disease free equilibrium is always globally asymptotically stable when R0, and 3. according to the choice of probability functions in the latent and temporary immune periods, establish that the disease always persists when R0 < 1 and an endemic equilibrium exists with different stability properties. In particular, the endemic steady state is at least locally asymptotically stable if the probability distribution in the temporary immunity is a decreasing exponential function when the duration of the latency stage is fixed or exponentially decreasing. It may become oscillatory under certain conditions when there exists a constant delay in the temporary immunity period. Numerical simulations are given to verify the theoretical predictions.

Bifurcations in a white-blood-cell production model.

We study the dynamics of a model of white-blood-cell (WBC) production. The model consists of two compartmental differential equations with two discrete delays. We show that from normal to pathological parameter values, the system undergoes supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles. We characterize the steady states of the system and perform a bifurcation analysis. Our results indicate that an increase in apoptosis rate of either hematopoietic stem cells or WBC precursors induces a Hopf bifurcation and an oscillatory regime takes place. These oscillations are seen in some hematological diseases.

Oscillations in cyclical neutropenia: new evidence based on mathematical modeling.

We present a dynamical model of the production and regulation of circulating blood neutrophil number. This model is derived from physiologically relevant features of the hematopoietic system, and is analysed using both analytic and numerical methods. Supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles are shown to exist. We make the estimation of kinetic parameters for dogs and then apply the model to cyclical neutropenia (CN) in the grey collie, a rare disorder in which oscillations in all blood cell counts are found. We conclude that the major cause of the oscillations in CN is an increased rate of apoptosis of neutrophil precursors which leads to a destabilization of the hematopoietic stem cell compartment.

Spatial effects in modeling pharmacokinetics of rapid action drugs.

We develop a model to describe the time course of plasma concentration of neuromuscular blocking agents used as anesthetics during surgery. This model, which overcomes the limitations of the classical compartment models routinely used in pharmacokinetics, incorporates spatial effects due to heterogeneity in the circulation: it takes the form of a dispersion equation on a circular domain, with a time-dependent leakage term. This term is fitted to the functional form desired once first-stage transients have died out. Comparisons are made with clinical data by adjusting three parameters in the model.

Complex dynamics and multistability in a damped harmonic oscillator with delayed negative feedback.

A center manifold reduction and numerical calculations are used to demonstrate the presence of limit cycles, two-tori, and multistability in the damped harmonic oscillator with delayed negative feedback. This model is the prototype of a mechanical system operating with delayed feedback. Complex dynamics are thus seen to arise in very plausible and commonly occurring mechanical and neuromechanical feedback systems. (c) 1995 American Institute of Physics.

Dynamical disease: Identification, temporal aspects and treatment strategies of human illness.

Dynamical diseases are characterized by sudden changes in the qualitative dynamics of physiological processes, leading to abnormal dynamics and disease. Thus, there is a natural matching between the mathematical field of nonlinear dynamics and medicine. This paper summarizes advances in the study of dynamical disease with emphasis on a NATO Advanced Research Worshop held in Mont Tremblant, Quebec, Canada in February 1994. We describe the international effort currently underway to identify dynamical diseases and to study these diseases from a perspective of nonlinear dynamics. Linear and nonlinear time series analysis combined with analysis of bifurcations in dynamics are being used to help understand mechanisms of pathological rhythms and offer the promise for better diagnostic and therapeutic techniques. (c) 1995 American Institute of Physics.